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ABSTRACT 



K^ , Terzan 5 is a globular cluster-like stellar system in the Galactic Bulge which 

^ . has been recently found to harbor two stellar populatio ns with different iron con 



tent and probably different ages (JFerraro et al. 1 120091 ) . This discovery suggests 



that Terzan 5 may be the relic of a primordial building block which contributed 
to the formation of the Galactic Bulge. Here we present a re-determination of 
the structural parameters (center of gravity, density and surface brightness pro- 
files, total luminosity and mass) of Terzan 5, as obtained from the combination 
of high-resolution {ESO-MAB and HST ACS-WFC) and wide-field (ESO-WFl) 
observations. We find that Terzan 5 is significantly less concentrated and more 
massive than previously thought. Still it has the largest collision rate of any 
stellar aggregate in the Galaxy. We discuss the impact of these findings on the 
exceptional population of millisecond pulsars harbored in this stellar system. 
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Subject headings: Globular clusters: individual (Terzan 5); stars: evolution 



1. INTRODUCTION 

Terzan 5 is commonly catalogued as a globular cluster (GC) located in the inner Bulge 
of our Galaxy. It is difficult to observe because it is heavily reddened, with an average color 
excess E{B -V) = 2.38 (Barbuy et al. 1998; Valenti et al. 2007). Not only is the reddening 
large, but it strongly depends on the line of sight (differential reddening; see Ortolani et al. 
1996, Valenti et al. 2007). Terzan 5 has an exceptionally large population of millisecond 
pulsars (MSPs). Indeed, the 33 MSPs detected so-far in Terzan 5 amount to about 25% of 
the entire sample of known MSPs in Galactic GCs (Ransom et al. 2005; see the updated list 



at http : //www . naic . edu/pf reire/GCpsr . html ) 



As part of a project (Ferraro et al. 2001, 2003; Cocozza et al. 2008) aimed at studying 
the properties of stellar populations harboring MSPs, we obtained a set of high-resolution 
images of Terzan 5 i n the K and J bands using the multi-conjugate adaptive optics (AO) 



demonstrator MAD ( iMarchetti et al. 1120071 ) temporally installed at the European Southern 
Observatory (ESO) Very Large Telescope (VLT). The {K, J — K) color- magnitude diagram 
(CMD) obtained from these observations led to the discovery of two well-defined red horizon- 
tal branch (HB) clumps, clearly separated in luminosity {6K ~ 0.3) and color (see Ferraro et 
al. 2009, hereafter F09; also see Figured]). A prompt spectroscopic follow-up demonstrated 
that the two populations have the same radial velocity (hence they belong to the same stel- 
lar systems) and their metal content is different: [Fe/H] ~ —0.2 and [Fe/H] ~ +0.3 for the 
fainter and the brighter group, respectively. These findings and the comparison with theoret- 
ical stellar isochrones con&m the existence of two distinct stellar populations in Terzan 5 and 
suggest that they possibly have been generated by two bursts of star formation with a time 
separation of a few (~ 6) Gyrs. While the age gap can be reduced by invoking a difference 



in the helium content of the two populations ( iD'Antona et al.ll2010[ ). the iron enrichment 
and the spatial segregation of the brightest clump, together with the extraordinary amount 
of MSPs found in Terzan 5, indicate that this system probably experienced a particularly 
troubled formation and evolutionary history (see Sect. 5). Terzan 5 is the ffist a GC-like 
system in the Galactic Bulge found to have a spread in the iron content and it could be the 



^ Based on observations collected at the European Southern Observatory, Cerro Paranal (during the 
Science Demonstration of MAD) and La Silla (under proposal 278.D-5067), Chile. Also based on observations 
with the NASA/ESA HST, obtained at the Space Telescope Science Institute, which is operated by AURA, 
Inc., under NASA contract NAS5-26555. 
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relic of one of the building blocks that contributed to the formation of the Bulge. Indeed the 
discovery might represent the observational evidence that even the innermost part of galactic 
spheroids form (at least partially) by t he accretion /rnergin g of small, previously formed and 



internally evolved stellar systems (e.g.. Ilmmeli et al.ll2004j ) 



In this paper we present the accurate re-determination of Terzan 5 structural parameters 
(surface density and surface brightness profiles, total luminosity, collision rate, etc.), obtained 
from a combination of high-resolution and wide-field observational data. These parameters 
provide basic information for a deeper understanding of the origin and the evolution of this 
puzzling system. 



2. THE DATA 

The photometric dataset used in the present work consists of high-resolution and wide- 
field images obtained, respectively, with MAD at the ESO-VLT and the Advanced Camera 
for Survey (ACS) on board the Hubble Space Telescope (HST), and with the Wide Field 
Imager (WFI) at the 2.2m ESO-MPI telescope. 

1. The MAD- dataset consists of a set of short (2-minute long), AO-corrected exposures 
secured through the K and J filters in August 2008, as part of a MAD science demonstration 
project (P.I. Ferraro). From all these images, we chose and analysed the highest quality 
ones: the full width at half maximum (FWHM) measured for the selected K and J images is 
0.1" and 0.24", respectively, just slighly larger than the diffraction limit. More importantly 
the FWHM is extremely stable over the entire (1' x 1') MAD field of view (FOV), fully 
demonstrating the potentiality of the multi-conjugate AO correction. 

2. The HST-dataset consists of deep ACS Wide Field Camera (WFC) images obtained 
through filters F606W (a broad V) and F814W (J), with total exposure times of 340 s and 
360 s, respectively (Prop. 9799, P.I. Rich). 

3. The Wide field-dataset consists of multiple V and I images (4 x 120 s exposures 
each) obtained with the ESO-WFI at La Silla (Chile) and retrieved from the ESO Science 
Archive (Prop. 278.D-5067(A); P.I. Bassa). The WFI is a mosaic of 8-CCD chips which 
combines wide-field (FOV =33' x 34') and reasonably high-resolution capabilities (pixel size 
of ~ 0.24"). The core of the cluster is roughly centered on CCD 7 and the images allow to 
sample the entire cluster extension (see Sect. H]). 

The point spread function (PSF) for each image has bee n modelled on several bright 



and isolated stars, by using the DAOPHOTH/PSF routine fIStetson 1 119871 ). Then PSF- 
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fitting photometry of the MAD dataset has been perfomed independently on the K and 
J best images, using DAOPHOTII/ALLSTAR. We then used DA0PH0TII/M0NTAGE2 
procedure to produce a stacked master frame combining the optical images in all the available 
filters (F606W and F814W for the ACS observations, and V and / for the WFI dataset). 
This method allowed us to obtain, for each optical dataset, a single high S/N image, cleaned 
from specific detector defects (i.e. hot pixels) and other spurious sources like cosmic rays 
etc. A master star list has been searched at a 4o" detection limit on this reference frame 



and it has then been used as input for ALLFRAME (jStetson Ill994l ). which simultaneously 
determines the brightness of the stars in all the frames, while enforcing one set of centroids 
and one transformation between all the images. Finally, the magnitudes obtained for each 
star have been normalized to a reference frame and averaged together, and the photometric 
error was derived from the standard deviation of the repeated measures. 

The star positions in the WFI sample were placed on the absolute astrometric system 
by using more than twenty thousand stars in common with the new astrometric Two Micron 
All Sky Survey (2MASS) catalogue!. Then the ~ 6000 stars in common between the WFI 
and the ACS datasets have been used as secondary astrometric standards for placing the 
ACS sample on the absolute astrometric system, and the same has been done for the ~ 
9000 stars in com mon between the AC S and the MAD datasets, following the procedure 



described, e.g., in iLanzoni et al. I ( 120071 ). The final astrometric accuracy of all the samples 
is of the order of ~ 0.2" both in right ascension (a) and in declination (5). The near-infrared 
instrumental magnitudes have been re ported to the 2MASS photometric system by using 



the stars in common with the catalog of IValenti et al. I (120071 ). The optical magnitudes have 



been transformed and homogenized by using a sample of stars in common with Ortolani et 
al. (1996). 

The CMDs obtained for the three datasets are shown in Figured! The two HB clumps 
discussed by F09 are clearly visible in the MAD dataset and partially distinguishable in the 
ACS sample, while the field contamination (mainly from the Galactic Bulge and Disk) is 
progressively more important at increasing distance (r) from the cluster center. The optical 
CMD clearly shows how Terzan 5 is strongly affected by differential reddening. We are 
currently constructing a detailed differential reddening map in the direction of this stellar 
system (Mucciarelli et al. 2010, in preparation), and preliminary results suggest that the 
reddening variation can be as high as 6 E[B — V) = 0.6. 



^Available at http://irsa.ipac.caltech.ed.u 
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3. CENTER OF GRAVITY 

F09 determined the center of gravity (Cgrav) of Terzan 5 by using the absolute positions 
of individual stars detected in the MAD sample. As a check, here we have recomputed Cgrav 
by exploiting the ACS dataset, with a cut in magnitude J = 21 needed to avoid spurious 
effects due to incompleteness in the very inner regions of the cluster. We found that, within 
the uncertainty (~ 0.5" in both a and 5), the two determinations are coincident, and we 
therefore confirm that gravity center of Terzan 5 is located at (F09): aj2ooo = 17*^ 48™ 4.85*^, 
^.72nnn = — 24 °46M4.6". This is ~ 1" north-we st from the photom etric centre quoted by 



Harris I ( 1l996l ). and ~ 1" north-east from that of ICohn et al. I ( 12002| ). 



To further check the reliability of o ur determination we also used a different approach 



This relies on the algorithm defined by ICasertano &: HutI ( 119851 ) to calculate the " density 



center" of a stellar system (for more details see Miocchi et al. 2010, in preparation). As 
a first step, for all the N stars located within a circle of a given radius, centered at an 
initially guessed position, the local (surface) density pi around the position of the ith star is 
evaluated as the inverse of its squared distance from the sixth nearest star. Then the cluster 
center is computed as the density-weighted average of the N star positions. The evaluation 
is repeated by centering the circle at the last found points and the iteration stops when 
the distance between the new and the previous determinaton is smaller than a given value 
(usually 0.01"). The density center thus computed agrees with the value quoted above well 
within the uncertainties. 

Finally, we also checked whether the two HB populations share the same gravity centre. 
To this purpose we have selected them from the ACS dataset, on the basis of the star 
position in the (J, V — I^ CMD (see Supplementary Figure 1 in F09). The projected spatial 
distribution of the two selected samples is plotted in Figure [2] and clearly shows that the 
bright HB population is more centrally concentrated than the faint one, in agreement with 
what found by F09 from the MAD dataset. The barycentres of the two HB populations seem 
to be different, with the faint HB centre being located ~ 3" south-east from that of the bright 
HB stars, which almost coincides with the cluster gravity centre quoted above. We stress 
however that, while the optical selection allowed us to improve the statistics and increase the 
radial coverage of the samples at most, it does not guarantee a proper separation of the two 
groups of stars. Hence, before confirming such a finding it is necessary to perform a more 
robust and clean selection of the two HB populations, based on (presently not available) 
near-infrared data covering a much larger area than the MAD FOV. 
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4. PROJECTED DENSITY AND SURFACE BRIGHTNESS PROFILES 

We have determined the projected density profile of Terzan 5 using direct star counts on 
the available datasets, which cover the entire radial extent of the cluster. For the innermost 
part of the profile we have exploited the high-resolution (MAD and ACS) datasets, while for 
r > 100" we have used the WFI sample and complementary data from the 2MASS survey, 
thus covering radial distances out to r = 1700". To avoid incompleteness biases, different 
limiting magnitudes have been adopted for the four datasets: K = 13, I = 20, / = 19, and 
K = 12.5 for the MAD, ACS, WFI, and 2MASS samples, respectively. By also adopting 
the colour cuts {V — I) > 3.4 and (J — K) > 1.3 we have excluded from the analysis the 
contribution of the stars belonging to the Galactic disk main sequence. With these limits 
we have computed the four portions of the density profile corresponding to each datasetjj 
In total more than 50,000 stars were used to c onstruct the star density profile. Following 



the procedure described in iFerraro et al.l ( 119991 ) we have divided the samples in concentric 
annuli centered on Cgrav? and each annulus has been split into an adequate number of sub- 
sectors. The number of stars lying within each sub-sector was counted, and the star surface 
density was obtained by dividing these values by the corresponding sub-sector areas. The 
stellar density in each annulus was then obtained as the average of the sub-sector densities, 
and the standard deviation was adopted as the uncertainty. Then, the radial annuli in 
common between the different samples have been used to shift and join the various portions 
of the profile. The overall projected density profile thus obtained is shown in Figure^] (empty 
squares), with the abscissae corresponding to the mid-point of each radial bin. The outermost 
(r ^ 175") measures from the WFI and 2MASS samples have an almost constant value, and 
their average has been used to estimate the Galactic Bulge and Disk contamination level. 
The subtraction of this background yields the profile shown in the figure as filled dots. The 
derived dens ity profile is well fit all over its radial extension by an isotropic, single-mass 



King model (JKing 1119661 ) with core radius re = 9.0", half-mass radius Vh = 31", tidal radius 
radius rt = 277" = 4.6', and intermediate concentration (c = 1.49). While the size of the 
core radius is consistent with the most recent determination (r^ = 7.9"; Cohn et al. 2002), 
the concentration is significantly smaller than that (c ~ 2) suggested by those authors, and 
the ratio between the core and the half-mass radius is a factor of two larger in our case. 

Exploiting the exceptional quality of the available datasets, we have also computed the 



■^While different limits have been (necessarily) adopted for the four datasets, the mass of the sampled 
stars is roughly constant (comparable to the that of the main sequence turn off stars). Hence any mass 
segregation effect on the four different portions of the profile is expected to be negligible. This is further 
confirmed by the agreeement found between the number density and the surface brightness profiles (see 
below). 
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surface brightness (SB) profiles by aperture photometry on the MAD and ACS images. The 
SB values were computed as the sum of the photon counts in each pixel, divided by the 
sampled area in any given radial annulus. The counts have been converted to a magnitude 
scale and then calibrated using a relation derived by performing aperture photometry on a 
number of high S/N isolated stars. The resulting SB profile, obtained for the inner ~ 100" 
from the center, after proper subtraction of the background level and of a few artefacts due 
to saturated stars, is shown in three different filters in Figure HI These profiles are well fit 
by the same King model derived from the projected density distribution, thus confirming 
that the samples of resolved stars used above are not affected by radial variations of the 
completeness and are properly selected. The values of the central SB measured in the three 
photometric bands are listed in Table 1, together with all the relevant parameters derived 
for Terzan 5. 



5. DISCUSSION 

The star density and SB profiles can be used to derive the integrated luminosity of 
the cluster. From the best-fit King model we estimate that the percentage of cluster light 
within regions of radius r = 15", 18", and 20" are roughly 30%, 36% and 40%, respectively. 
Using aperture photometry on the MAD images, we obtain integrated-light values of K{r < 
15") = 3.44, K{r < 18") = 3.3 and K{r < 20") = 3.2 mag, respectively. Adopting the 
color excess E{B — V^) = 2.38, the distance modulus (m — M)o = 13.87 (Valenti et al. 
2007, corresponding to a distance oi d = 5.9 ± 0.5 kpc) and the bolometric correction 
BCk = 2.4 appropriate for a population of intrinsic color (J — K)q = 0.8 (see Montegriffo et 
al. 1998), we estimate that the corresponding bolometric luminosity in the considered regions 
is: Lboi(r < 15") = 3 x IO^Lq, Lboi(r < 18") = 3.4 x IO^Lq and Lboi(r < 20") = 3.7 x IO^Lq. 
Considering the fraction of light sampled in each region we find that the total luminosity of 
the system is Lboi = 9.5 ± 0.3 x IO^Lq. 

An independent estimate of the total luminosity of the stellar system can be derived 



from its stellar population, by using a simple relation (JRenzini fc Buzzoni Ill986l ) linking the 
number of stars (Nj) observed in a given post-main sequence evolutionary stage j and the 
luminosity of the entire parent cluster (Lt)'- 

Nj = B X tj X Lt, (1) 

where B is the specific evolutionary flux (for intermediate/old stellar populations B = 2 x 
10~^^ stars yr~^LQ^) and tj is the duration of the evolutionary stage. The number of HB 
stars counted in the two clumps by F09 is quite large: a total of about 1300 (with 800 



and 500 belonging to the faint and the bright HB clumps, respectively). This population 
is compa rable to (or even larg er than) that o bserved in the largest Ga lactic GCs, like 47 
Tucanae (JBeccari et al. Il2006al ) and NGC 6388 (JDalessandro et al. II2008I ). and suggests that 
the overall size of Terzan 5 (in terms of luminosity and mass) is comparable to that of these 
systems. For a quantitative estimate, we insert the observed number of HB stars in the above 
relation and adopt tnB = 10^ yr. This provides a luminosity of 4 x IO^Lq and 2.5 x IO^Lq 
for the two parent populations, and a total luminosity of 6.5 x IO^Lq for the entire stellar 
system. This estimate, which is distance and reddening independent, is quite consistent 
with the previous one, thus confirming that Terzan 5 has a considerable total luminosity 
(hereafter we adopt the average value Lboi = 8 x 10^ Lq), significantly higher than previously 
thought. By comparison, adopting the values of distance and reddening quoted above and 
a bolometric correction Lboi — 1.4 Ly, the tot al bolorn e tric lu minosity corresponding to the 
integrated magnitude (Vt = 13.85) quoted by [Harris I ( 119961 ) would be only Lboi — lO^L©. 
The discrepancy is most probably due to the strong (differential) reddening affecting the 
system, especially in the optical bands. This effect is greatly reduced for our new estimate, 
since it is based on the observed i^-band integrated r nagnitude and t he number of HB stars. 
By assuming a mass-to-light ratio M/Lboi = 3 (e.g., lMarastonlll998l ). the total stellar mass 
of this system is Mt ~ 2 x lO^M©. 



Verbunt fc Hut I (Il987l ) first suggested that the collision rate of Terzan 5 is the highest 



among the Galactic GCs. We ca n now re-compu t e this quantity by adopting the newly 
determined parameters. Following IVerbunt fc Hut I ( 119871 ). the coUisional parameter (F) for 
a King-model or virialized system can be computed as: F oc po ^ r^'^, where pn is t he cen tral 
mass density. By using the values obtained above and equation (7) of iDjorgovski I ( 1l993l ). we 
find that the collision parameter of Terzan 5 is between 5 and 10 times higher than that of 
Filler 1 and of other massive clusters for w hich structural parameters have been recently re 



determined (NGC6 388, NGC6266, 47Tuc; IPalessandro et al. I l2008l : iBeccari et al. I l2006b 



Mapelli et al.l l2006l respectively). Hence we confirm that, even with the new structural 
parameters (suggesting a lower concentration and a larger mass than previously thought), 
Terzan 5 still has the largest known collision rate of any stellar aggregate in the Galaxy. 

The co-existence of two stellar populations with different iron content (and probably 
ages) suggests that the original mass of Terzan 5 was significantly larger in the past than 
observed today, large enough to retain the iron-enriched gas that, otherwise, would have been 
ejected out from the system by the violent supernova (SN) explosions. Indeed, the smallest 
systems with solid evidences of a spread in the iron content (and ages) are significantly more 
massive than GCs: the dwarf spheroidal satellites of the Milky Way typically have masses of 
~ IO'^Mq (Strigari et al. 2008; see also iBattaglia et al.ll2008l ) and, following recent chemo- 
dynamical models well reproducing the observations, their initial masses amounted to a few 
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10^ Mq (JRevaz et al.ll2009f). While a lower l i mit o f ~ lO^M© for the proto-Terzan 5 could 
also be hazarded following iBaumgardt et al.l (120081 ) . more detailed and extensive simulations 
are needed to firmly determine the smallest total mass necessary to retain the SN ejecta. 

The exceptionally high metallicity regime of the two stellar populations found in Terzan 5 
also suggests a quite efficient enrichment process, that could have a relevant role in the origin 
of its population o f MSPs. In particular , both the iron and the [a/Fe] abundance ratios mea- 
sured in Terzan 5 ( lOriglia fc Rich II2004J . Rich et al. 2010, in preparation) show a remarkable 
similarity with those of the Bulge stars. This strongly suggests that these two structures 
shared the same sta r formation and chemi c al enrichment processes. The many ob servations 



of Bulge stars (e.g., iMelendez et al. II2008I : lOrigha et al.l 120081 : iRvde et al. II2009I . and refer- 



ences therein) indicate that they are all characterized by an old age, a high (close to solar) 
average metallicity [Fe/H], and an [a/Fe] ratio which is enhanced (due to SNII enrichment) 
up to a metallicity [Fe/H]~ 0. These constraints suggest a scenario where the dominant 
stellar population of the Bulge formed early (thus explaining the old agcl), rapidly and with 
high efficiency (from a gas mainly enriched by SNII, thus explaining t he [a/Fe] enhance- 



ment up to high iron c ontentp. Also chemical evolution models (e.g., iBallero et al.l 120071 : 
Mc William et al.ll2008l ) indicate that the abundance patterns observed in the Bulge require 
a quite high star formation efficiency and an initial mass function flatter than that in the 
solar neighbourhood, thus to rapidly enrich the gas up to about solar metallicity through 
an exceptionally large amount of SNII explosions. The assumption of a similar scenario for 
Terzan 5 would naturally explain its extraordinary population of MSPs, since the expected 
high number of SNII would produce a large population of neutron stars, most of which would 
have been retained by the deep potential well of the massive proto-Terzan 5 system. Then 
the high collision rate could have favoured the formation of binary systems containing neu- 
tron stars and promoted the re-cychng process that flnally generated the large population of 
MSPs now observed in Terzan 5. If su ch a scenario is corre ct, many more MSPs still wait to 
be discovered in this system (see also JRansom et al.ll2005l ). the 33 known objects probably 
being just the tip of the iceberg. Future deeper pulsar searches of Terzan 5, perhaps with 
larger telescopes such as the Square Kilometer Array, will shed additional light on the nature 
of this system. 



This research was supported by the Agenzia Spaziale Italiana (under contract ASI- 



Additional episodes of sta r formatioii mainly confined in the innermost (~ 100 pc) region could eventually explain the 
presence of younger stars fe.g. jBlum et al.ll2003l : iFieer et al.ll2004r . 

The [Q:/Fe]-[Fe/H] relation shows a down-turn at a value of [Fe/H] which depends on the star formation rate; the higher 
the latter, the higher the metallicity at which the down-turn occurs. Such a value is [Fe/H]~ —1 in the Old Halo/Disk, while 
it is significantly higher ([Fe/H]~ 0) in the Bulge, testifying a much higher star formation rate in this dense environment. 
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Fig. 1. — Color-magnitude diagrams for the three datasets used in the paper and for different 
radial regions of the cluster (see labels). The double HB discovered by F09 is clearly visible 
in the MAD near-infrared CMD, while it is more difficult to distinguish in the ACS optical 
plane (but see Supplementary Figure 1 in F09). The WFI sample is dominated by field star 
contamination. 
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Fig. 2. — Map of the bright {left panel) and faint {right panel) HB populations selected 
from the ACS dataset on the basis of the star position in the (/, V — I) CMD. The star 
coordinates {a and 6) are plotted with respect to those of the cluster gravity centre {ac and 

5r). 
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Fig. 3. — Newly determined star density profile of Terzan 5, obtained from resolved star 
counts in the combined photometric dataset. Empty squares represent the observed profile, 
while solid dots are obtained after subtraction of the field background density (marked with 
the dotted line). The best-fit single-mass King model is shown as a solid line and the 
corresponding structural parameters (core radius Vc, concentration c and tidal radius r^) are 
labelled in the figure. The lower panel shows the residuals between the observations and the 
fitted profile at each radial coordinate. 
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Fig. 4. — Surface brightness profile of Terzan 5 in three different filters, after background 
subtraction: K band (upper panel), obtained from the MAD dataset, / and V bands (middle 
and lower panels), from the ACS sample. In each panel the solid line corresponds to the 
best- fit King model shown in Fig.O 
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Table 1. Structural parameters for Terzan 5 



Center of gravity 

Reddening^ 
Distance^ 
Core radius 
Concentration 
Total luminosity 
Total mass 
Central mass density 
Central i^-band SB 
Central /-band SB 
Central l^-band SB 



aj2ooo = 17M8-4.85« 

5j2ooo = -24°46'44.6" 

E{B-V) = 2.38 ±0.055 

rf = 5.9±0.5kpc 

r^ = 9" = 0.26 pc 

c= 1.49 

Lboi ^8x10^^0 

Mt ~ 2 X IO^Mq 

po ^4.1 X lO^Mo/pc^ 

fJ'K{0) = 9-85 mag/arcsec^ 

/i/(0) = 15.87 mag/arcsec^ 

/iv'(O) = 20.54 mag/arcsec^ 



Note. — tfrom IValentLetalJ (120071 ) 



